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I.  INTRODUCTION 


1.1  BACKGROUND 

The  X-wing  concept  is  a  novel  approach  to  the  problem  of 
combining  VTOL  and  hover  capability  with  efficient  high-speed 
cruise  operation  ll].  The  concept  uses  a  rotating  X-platform 
wing  to  operate  like  a  conventional  helicopter  for  VTOL,  hover 
and  low-to-moderate  speed  operation;  the  rotor  may  also  be 
stopped  when  at  sufficient  forward  speed,  allowing  the  aircraft 
to  cruise  as  a  fixed-wing  vehicle. 

A  key  technology  element  of  this  concept,  as  developed  by 
the  Lockheed  California  Aircraft  Company  in  their  preliminary 
design  study  L2],  is  the  use  of  a  very  rigid,  circulation- 
controlled  rotor.  Lift  and  control  moments  are  maintained  in 
both  rotary-  and  fixed-wing  flight  by  ejection  of  low-energy  air 
from  the  trailing  edge  of  the  blade. 

— *  The  original  purpose  of  this  contract  was  to  prepare  a 
flight  test  plan  for  the  proposed  X-wing  demonstrator  using 
system  identification  to  extract  useful  math  models  from  the 
flight  test  data.  Since  the  original  statement  of  work  was 
submitted,  however,  the  scope  of  this  study  changed.  An  X-wing 
program  decision  has  been  reached  not  to  carry  the  Lockheed 
X-wing  configuration  to  flight  test.  Therefore,  this  task  was 
modified  to  provide  a  demonstration  of  the  feasibility  of  using 
system  identification  techniques  to  extract  low-order  math  models 
from  time  history  data  from  a  detailed  X-wing  rotor  simulation 
CREXORJ.  -  - 

The  X-wing  concept  has  been  simulated  by  Lockheed  using  a 
general-purpose  helicopter  simulation  program  known  as  REXOR 
13] .  REXOR  models  the  rotor  in  a  rotating  frame  of  reference  and 
models  each  biade  separately.  This  program  has  been  installed  by 


SCI  (Vt)  on  the  CDC-7600  computer  system  at  NASA-Ames  Research 
Center,  Moffett  Field,  CA,  as  a  basis  for  this  study. 

1.2  METHOD  OF  APPROACH 

The  REXOR  simulation  program,  as  implemented  on  the  CDC  7600 
at  NASA-Ames,  provided  the  data  base  from  which  the  X-wing  math 
models  were  to  be  identified.  Numerous  simulation  runs  were 
made,  and  the  data  files  produced  were  written  to  tape  and 
transported  to  SCI  (Vt)  for  anaysis  on  SCI's  VAX  computer 
system.  Data  used  represented  hover  and  100-foot-per-second 
flight  conditions  for  collective,  longitudinal  and  lateral 
cyclic,  and  tail  rotor  collective  inputs.  A  reduced-order  math 
model  of  the  X-wing  was  then  identified  using  SCI's  Optimal 
Subset  Regression  (OSR)  program.  The  simplified  model  describes 
the  fuselage  degrees  of  freedom  and  the  rotor  hub  moments. 

One  of  the  problems  anticipated  in  identifying  a  low-order 
model  from  he  REXOR  data  is  that  the  REXOR  model  includes 
high-frequency  (u>^  >  fiR)  aeroelastic  modes  which  are  of 
little  value  to  the  flight  controls  designer.  For  this  reason, 
it  was  expected  that  prefiltering  the  data  to  remove  this 
unnecessary  information  would  result  in  better  identification  of 
the  simplified  model,  since  this  reduced-order  model  cannot 
account  for  the  high-frequency  effects  seen  in  the  data. 

1.3  REPORT  SUMMARY 

The  second  section  of  this  report  presents  a  description  of 
the  X-wing  configuration.  The  intent  of  this  description  is  to 
provide  background  information  for  those  features  of  the  rotor 
design,  the  vehicle's  dynamic  characteristics,  and  the  flight 
control  system  which  have  an  impact  on  the  modeling  approach. 
Section  III  describes  the  methodology  for  generating  simulation 
time  history  data  with  REXOR,  the  X-wing  nonlinear,  high-fidelity 
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simulation  math  model.  Section  IV  describes  the  developed 
methodology  for  extracting  reduced-order  models  of  the  X-wing 
configuration  from  REXOR  simulation  data.  Report  conclusions  and 
recommendations  are  presented  in  Section  V. 


II.  X-WING  DESCRIPTION 


This  chapter  presents  a  description  of  the  Lockheed  X-wing 
configuration.  The  intent  of  this  description  is  to  provide 
background  information  for  those  features  of  the  rotor  design, 
the  vehicle's  dynamic  characteristics,  and  the  flight  control 
system  which  have  an  impact  on  the  modeling  approach.  The  X-wing 
configuration  is  shown  in  Figure  2.1,  and  its  geometric 
parameters  are  listed  in  Table  2.1. 

2.1  ROTOR  DESIGN 

The  principal  feature  of  the  X-wing  aircraft  is  a  single, 
stoppable  rotor  on  top  of  the  fuselage,  which  enables  the 
airplane  to  operate  like  a  helicopter  with  the  rotor  turning,  or 
like  a  fixed-wing  aircraft  with  the  rotor  stopped  in  any 
X-position  relative  to  the  airplane's  centerline.  The  only  other 
lifting  surfaces  on  the  aircraft,  as  presently  conceived,  are  on 
the  empennage.  In  the  helicopter  mode,  forward  propulsion  is 
provided  by  both  rotor  axis  tilt  and  auxiliary  engines.  In  the 
stopped-rotor  mode,  propulsion  is  provided  directly  by  the 
engines.  The  rotor  is  a  key  advancement  in  the  state  of  the  art; 
it  is  a  very  rigid,  hingeless  rotor,  much  stiffer  than  the 
so-called  "rigid"  rotors  used  on  some  helicopters,  and  able  to 
support  the  vehicle  and  provide  maneuvering  capability  in 
stopped-rotor  cruise.  Even  the  requirement  for  a  collective 
hinge  has  been  eliminated  by  use  of  circulation  control  (CC)  to 
regulate  blade  lift.  The  CC  concept  involves  a  stream  of  high 
pressure  air  exhausted  from  the  trailing  edge  of  the  elliptical 
blade  airfoil  section  which  fixes  the  rear  stagnation  point 
(through  the  Coanda  effect),  and  results  in  a  net  circulation 
about  the  airfoil.  Hence,  the  common  name  for  the  system  is 
circulation-controlled  rotor  (CCR).  The  X-wing  uses  the  CC 
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Table  2.1 
X-W ing  Geometry 


PARAMETER 

DESCRIPTION 

VALUE 

UNITS 

*F 

c.g.  location  {fuselage  axes) 

0.31 

ft 

Hr 

Height  of  hub  above  fuselage 
ref. 

1.13 

ft 

*TR 

Distance  aft  to  tail  rotor 

14.80 

ft 

Ns 

Distance  aft  to  vert,  stabilizor 

15.80 

ft 

Ns 

Oi stance  aft  to  horiz.  stabilizor 

19.10 

ft 

Hvs 

Height  of  vert,  stabilizor 

6.88 

ft 

htr 

Height  of  tail  rotor 

3.67 

ft 

Zp 

Dist.  of  thrustor  below 
fuselage  ref. 

-  0.58 

ft 

r*R 

Radius  of  main  rotor 

12.50 

ft 

/ 


concept  in  a  mode  wherein  high  advance  ratio  capability  is 
achieved  by  using  dual-slotted  airfoils  for  useful  lift 
production  in  the  reverse  flow  regime.  This  mode  is  identified 
as  the  "reverse  blowing  circulation-controlled  rotor  (RBCCR).” 

Theoretical  and  experimental  studies,  many  funded  by  ONR  and 
NASC,  have  confirmed  the  feasibility  of  the  RBCCR  concept. 
Aircraft  characteristics  are  satisfactory  in  both  performance  and 
control,  according  to  analysis  and  small-scale  wind  tunnel  data 
11].  Aerodynamic  characteristics  of  the  aircraft  with  the  rotor 
turning  and  with  it  locked  have  been  predicted  using  a  number  of 
Lockheed  and  government  programs  [2,3J.  The  characteristics 
combine  those  of  helicopters,  fixed-wing,  and  powered-lift 
aircraft,  and  are  therefore  dependent  on  angle  of  attack,  blowing 
pressure,  and  rotor  RPM.  Analyses  of  thrust  (lift)  and  moment, 
both  critical  to  performance  and  longitudinal  control,  show 
adequate  trim  capability  and  acceptable  static  margin  with  the 
addition  of  a  small  (6  ft^)  horizontal  tail.  Of  special 
interest  are  neutral  point  and  trim  shifts  at  transition.  The 
ability  to  use  the  fixed-wing  for  both  lift  and  pitch  control 
introduces  a  "direct  lift  control"  capability  and  reduces  the 
role  of  the  horizontal  tail  to  that  of  minimizing  trim  drag. 

2.2  DYNAMIC  CHARACTERISTICS 

Flight  dynamics  studies  have  been  performed  L  3  ]  to  analyze 
the  aircraft's  dynamic  modes  of  motion,  both  rigid-body  and 
aeroelastic,  in  various  flight  conditions.  The  predominant 
dynamic  effect,  in  the  helicopter  mode,  is  the  inertial 
cross -coupling ,  which  causes  unacceptable  handling  qualities  if 
not  compensated  by  an  automatic  control  system.  A  principal  area 
of  study  in  X-mode  operation  (rotor-stopped)  is  static  stability 
and  control  authority.  The  conversion  process  contains  its  own 
areas  of  concern,  particularly  oscillatory  loads  caused  by  a 
slowly  turned  rotor  and  their  effect  on  dynamics  and  handling 
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qualities.  An  RPM  load-factor  schedule  must  be  defined  for  the 
conversion  and  reconversion  maneuvers. 

Each  of  the  three  major  flight  regimes  (helicopters, 
X-cruise,  and  transition)  has  been  analyzed  using  linear  and 
nonlinear  models,  and  no  fundamental  flaws  in  the  concept  have 
been  found.  The  rotor  aeroelastic  modes  are  of  very  high 
frequency  (100  rad/sec)  due  to  the  extreme  rigidity  of  the  rotor, 
and  do  not  couple  with  the  aircraft  rigid-body  modes.  The 
eigenanalysis  discussed  in  Ref.  3  predicts  a  static  instability 
in  the  helicopter  mode.  It  appears  that  a  full-time, 
f ai 1 -operatonal  augmentation  system  is  required  for  stability 
reasons,  as  well  as  for  control  of  rotor  cross -coupl i ng .  The 
approximately  1.2  sec  time  to  double  pitch  amplitude  discussed  in 
Ref.  3  is  not  acceptable  under  such  handling  qualities  criteria 
as  MIL-F-83300. 

2.3  FLIGHT  CONTROL  SYSTEM 

Control  of  the  X-wing  configuration  is  based  on  the 
regulation  of  the  air  pressure  level  in  the  circulation  control 
(CC)  plenums,  a  lift  fan  in  the  vertical  tail  and  aerodynamic 
empennage  controls.  The  CC  pneumatic  blowing  for  each  blade  is 
regulated  as  a  function  of  blade  position  to  achieve  collective 
and  cyclic  control. 

Two  types  of  augmentation  systems  have  been  considered  for 
te  aircraft,  referred  to  as  low-gain  and  high-gain  systems.  Both 
are  aircraft  rate  command  systems,  but  the  high-gain  system 
features  a  high-gain  rotor  moment  feedback  loop  within  the 
aircraft  rate  loop  that  drives  the  rotor  controls  to  produce  a 
commanded  moment  and  null  out  disturbances  from  turbulence  or 
inertial  cross -coupl ing .  The  X-wing  aircraft,  as  modeled  for 
this  study,  uses  the  high-gain  control  system  for  the  pitch  and 
roll  degrees  of  freedom  (i.e.  cyclic  augmentation). 


9 


The  cyclic  control  augmentation  system  is  shown  in  Figure 
2.2.  Rotor  hub  moments  (MH  and  L^)  are  fed  back  at  high  gain 
to  remove  external  aerodynamic  effects  on  the  rotor.  Phasing 
networks  were  included  in  the  initial  design  for  the  moment 
feedback  loop,  but  analysis  showed  negligible  benefit  so  the 
feedback  phase  angle  was  set  to  zero  in  the  REXOR  simulation. 
Pitch  and  roll  rates  are  also  fed  back  to  augment  damping. 

Normal  acceleration  (Az)  feedback  was  not  implemented. 

In  the  forward  path,  control  input  phasing  is  used  to 
compensate  for  the  phase  shift  in  rotor  moment  with  advance 
ratio.  This  phase  angle  must  be  adjusted  for  rotor  speed,  and  is 
equal  to  60°  for  maximum  RPM.  Since  all  runs  for  this  study 
were  made  at  100%  RPM,  the  phase  angle  was  a  constant.  Forward 
loop  integration  is  used  to  stabilize  the  aircraft,  and  the 
acutators  control  the  area  of  cosine  and  sine  blowing  values, 

A^b  and  B^g.  This  gives  a  total  blowing  pressure  of 

PB  =  2116.  x  RpQ  (1.  -  A1B  cosi -  B1B  sin^R)  psf 

where 

Rpo  =  collective  pressure  ratio. 

The  gains  and  other  parameters  used  in  the  cyclic  control 
system  are  listed  in  Table  2.2. 
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Table  2.2 

Flight  Control  System  Parameters 


PARAMETER 

DESCRIPTION 

VALUE 

UNITS 

Kp 

Longitudinal  stick  gain 

1.0 

ft/ft 

Kr 

Lateral  stick  gain 

1.0 

ft/ft 

Kq 

Pitch  rate  feedback  gain 

2.52 

ft-sec/rad 

wq 

Pitch  rate  filter  freq. 

120. 

rad/sec 

q 

Pitch  rate  filter  damp. 

0.7 

— 

Kp 

Roll  rate  feedback  gain 

2.8 

ft-sec/rad 

Wp 

Roll  rate  filter  freq. 

120. 

rad/sec 

P 

Roll  rate  filter  damp. 

0.7 

— 

K 

Pitch  moment  feedback  gain 

2.1  x  10-4 

ft/ in- lb 

W 

Pitch  moment  filter  freq. 

120. 

rad/sec 

kl 

Roll  moment  feedback  gain 

3.2  x  10-4 

ft/in-lb 

WL 

Roll  moment  filter  freq. 

120. 

rad/sec 

Kir 

Cosine  blowing  gain 

0.7 

(ft)-l 

kir 

Sine  blowing  gain 

0.7 

(ft)-l 

T  c 

Actuator  time  constant 

0.033 

sec 

^  c 

Rotor  control  input  phase 

1.047 

rad 
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III.  GENERATION  OF  SIMULATION  DATA 


3.1  OVERVIEW 

The  first  task  in  this  study  was  to  produce  time  history 
data  representative  of  the  X-wing  design  proposed  by  Lockheed. 
This  was  done  using  the  REXOR  simulation  program  developed  by 
Lockheed.  The  simulation  was  implemented  and  run  on  the  CDC  7600 
computer  system  at  NASA-Ames  for  this  study. 

3.2  REXOR  DESCRIPTION 

The  REXOR  (Revised  and  Extended  Roto_r  Senior)  program  is  a 
general-purpose  simulation  developed  for  analysis  of  rotorcraft. 
The  version  of  REXOR  implemented  at  NASA-Ames  has  been  modified 
to  simulate  the  X-wing  configuration  designed  by  Lockheed,  and  a 
data  deck  representing  the  X-wing  was  also  available  at 
NASA-Ames.  The  basic  REXOR  program  is  documented  in  Ref.  3.  No 
documentat ion  exists  for  the  X-wing  version. 

The  standard  REXOR  program  was  written  to  simulate  a  single, 
two-  or  four-bladed,  hingeless -rotor  helicopter  with  a  gyro 
control  on  the  rotor.  The  program  also  has  the  ability  to 
simulate  conventional  controls.  The  X-wing  version  contains 
modifications  to  simulate  the  pneumatic  controls  used  in  that 
design,  as  well  as  simulating  the  transition  to  fixed-wing 
flight . 

The  REXOR  math  model  was  developed  to  provide  accurate 
analysis  of  rotor  characteristics.  This  has  been  achieved  using 
a  blade-element  approach,  in  which  the  blade  aerodynamics  and 
kinematics  are  computed  as  a  summation  of  finite  elements  of  each 
blade.  This  approach  is  more  complex  than  the  classical  tip-path 
plane  models  commonly  used,  but  it  allows  in-depth  harmonic 
analysis  of  the  aeroelastic  rotor  dynamics.  However,  it  is  also 


cumbersome  to  operate,  and  is  computationally  expensive.  As  an 
example  of  this,  one  2-sec  simulation  run  by  SCI  (Vt)  on  the  C DC 
CYBERNET  system  took  1030.417  sec  of  CPU  time.  Obviously,  this 
is  not  an  efficient  program  for  handl ing - qual i t  ies  and  flight 
control  design  studies  which  require  repeated  runs,  but  do  not 
require  extreme  accuracy  in  the  modeling  of  high-frequency 
aeroelastic  deflections  of  the  rotor.  By  identifying  a  lower 
order  model  of  the  X-wing  aircraft,  a  tool  can  be  developed  which 
would  allow  the  flight  control  designer  to  make  many  runs 
efficiently,  and  produce  the  data  needed  for  control  system 
design  and  handling  qualities  analysis. 

The  math  model  used  by  the  REXOR  program  also  complicates 
the  process  of  determining  a  low-order  model  in  that  it  is  not 
feasible  to  extract  stability  derivatives  by  the  perturbation 
method  commonly  used  with  fixed-wing  aircraft  simulations.  This 
is  because  the  aerodynamics  of  the  rotor  cannot  be  determined 
statically  from  a  blade-element  model.  An  entire  revolution  of 
the  rotor  must  be  computed  in  order  to  get  the  averaged  effects 
of  the  rotor  states.  For  this  reason,  system  identification  is 
the  logical  way  to  generate  reduced-order  models  from  the  REXOR 
math  model. 

3.3  TEST  DATA  GENERATION 

The  control  inputs  used  for  generating  the  time  histories  to 
be  identified  were  chosen  to  excite  any  modes  within  the 
frequency  range  of  interest  for  handling -qualities  analysis.  The 
inputs  to  all  four  pilot  controls  (main  rotor  collective,  tail 
rotor  collective,  longitudinal  cyclic  stick  and  lateral  cyclic 
stick)  consisted  of  a  series  of  four  doublets.  Each  doublet  had 
twice  the  period  and  half  the  magnitude  of  the  previous  one. 

These  inputs  are  shown  in  Figure  3.1.  Note  that  the  tail-rotor 
collective  input  was  the  same  as  the  main-rotor  collective 
blowing  input,  but  the  time  history  shown  was  the  actuator 
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Figure  3.1  Control  Inputs 


output.  The  actuator  for  collective  blowing  was  not  simulated  in 
the  code. 

After  analyzing  the  generated  REXOR  simulation  responses,  it 
was  apparent  that  very  little  pitch  and  roll  rate  response  was 
achieved  with  the  cyclic  stick  inputs.  This  is  shown  in  Figures 
3.2  and  3.3.  The  cyclic  data  were,  however,  included  for  system 
identification,  since  the  feedbacks  and  an  initial  miss  trim 
caused  the  cyclic  blowing  (A1B  and  B1B)  to  drive  the 
simulation  dynamics  throughout  these  maneuvers- 

The  miss  trim,  primarily  seen  as  a  hub  moment  transient,  is 
most  likely  due  to  an  initialization  problem  with  REXOR.  When 
REXOR  was  first  loaded  onto  the  NASA-Ames  CDC  7600  computer, 
r^merous  initialization  problems  ( i . e .  ,  undefined  variables)  were 
encountered.  The  approach  adopted  at  the  time  was  initially  to 
set  all  variables  to  zero.  This  across-the-board  initialization 
of  all  variables  to  zero  could  cause  problems  with  the  hub  moment 
feedback  filter.  The  filter  was  probably  initialized  correctly 
in  the  trim  mode,  but  then  the  filter  output  could  have  been 
zeroed  after  switching  the  simulation  to  the  run  mode.  From  the 
standpoint  of  identifying  unaugmented  airframe  math  models,  which 
is  the  objective  of  this  study,  this  miss  trim  situation  provides 
acceptable,  albeit  unnatural,  control  excitation.  In  conclusion, 
better  documentation  for  the  X-wing  version  of  REXOR  would  have 
helped  this  problem. 

The  test  cases  used  for  this  study  included  two  flight 
conditions:  hover  and  100-feet-per -second  forward  velocity. 

Both  cases  were  run  with  the  aircraft  in  the  helicopter  mode  of 
flight.  For  each  flight  condition,  four  runs  were  made,  one  for 
each  control  input.  Each  of  the  runs  made  had  a  duration  of  10 
seconds.  These  test  cases  are  summarized  in  Table  3.1. 

For  the  identification  processing,  the  four  runs  for  each 
flight  condition  were  concantenated .  This  allows  simultaneous 
identification  of  the  effects  of  all  four  control  inputs. 


TIME  (SECOND 


TIME  (SECONDS)  |  “  TIME  (SECONDS) 


Table  3.1 
Test  Cases 


RUN  I.D. 

FLIGHT  CONDITION 

CONTROL  INPUT 

1004.02 

Hover 

Collective  blowing 

1014.01 

Hover 

Longitudinal  stick 

1024.01 

Hover 

Lateral  stick 

1034.01 

Hover 

Tail  rotor  collective 

1104.02 

VT  =  100  FPS 

Collective  blowing 

1114.01 

VT  =  100  FPS 

Longitudinal  stick 

1124.01 

VT  =  100  FPS 

Lateral  stick 

1134.01 

VT  =  100  FPS 

Tail  rotor  collective 

Prior  to  making  production  data  runs  with  REXOR,  the  output 
subroutine  was  modified  to  store  the  data  reqired  for  the 
identification  process.  The  variables  which  were  output  and 
stored  are  listed  in  Table  3.2. 

3.4  POST -SIMULATION  DATA  PROCESSING 

Before  using  the  REXOR  simulation  data  with  the 
identification  program,  it  was  checked  for  consistency  between 
redundant  variables  (i.e.,  does  dQ/dt  =  Q).  The  REXOR 
simulation  data  were  also  filtered  to  remove  mult iple-per -rev 
dynamics  associated  with  the  structural  degrees  of  freedom. 

3.4.1  Kinematic  Consistency 

Examination  of  the  REXOR  output  data  showed  a  discrepancy 
between  the  time  histories  of  the  pitch  and  roll  rates  and  their 
derivatives.  Figures  3.4  through  3.6  show  the  comparison  between 
the  angular  accelerations  output  by  REXOR  (PDOT,  QDOT,  and  RDOT ) 
and  the  numerical  derivatives  of  the  angular  velocities  (dP/dT, 
dQ/dT,  and  dR/dT) . 

Further  examination  of  the  data  showed  that  the  angular 

velocities  P  and  Q  were  consistent  with  the  Euler  angle 

•  • 

rates,  4  and  Q  .  This  was  shown  to  be  the  case  by  comparing 
time  histories  of  P  and  Q  with  reconstructed  time  histories 
calculated  from  and  9  .  These  comparisons  are  shown  in 
Figure  3.7  and  3.8. 

Since  P  and  Q  were  consistent  with  4  and  0  ,  but  not 
with  P  and  Q,  the  numerical  derivatives  of  P  and  Q  were 
used  for  the  system  identification.  The  numerical  derivative  of 
R  was  also  used,  although  it  was  virtually  identical  to  the  R 
output  by  REXOR. 


Table  3.2 

REXOR  Output  Variables 


CHANNEL 

FORTRAN  NAME 

VARIABLE 

UNITS 

CHANNEL 

FORTRAN  NAME 

VARIABLE 

UNITS 

1 

ST1ME 

am 

sec 

21 

AUXE(4) 

8  IB 

non-dim 

2 

UFO 

ft/sec2 

22 

PS1R 

R 

rad 

3 

VFO 

&■ 

ft/sec2 

23 

(XCS-XCST)*12. 

xcs 

inches 

4 

HFD 

ft/sec2 

24 

(YCS-YCST)*12. 

Ycs 

inches 

5 

PRO 

p 

2 

rad/ sec 

25 

RPO 

8P0 

non-dim 

6 

quo 

Q 

? 

rad/sec 

26 

THOTR 

OTR 

rad 

7 

RRO 

R 

2 

rad/sec 

27 

REFS1G 

.1  *lr*R 

non-dim 

8 

PH  ID 

♦ 

rad/sec 

28 

XC1 

XC1 

inches 

9 

THETD 

6 

rad/sec 

29 

YC1 

YC1 

Inches 

10 

AUXE(l) 

A1B 

non-dim 

30 

V!2,l,l)*12. 

A21 

inches 

11 

AUXE(2) 

8ib 

non-dim 

31 

V(2,2,l)*12. 

A21 

In/sec 

12 

UF 

uf 

ft/sec 

32 

V(2,3.1)*12. 

A21 

in/sec2 

13 

VF 

vf 

ft/sec 

33 

V(2,l,2)*12. 

A22 

Inches 

14 

WF 

“f 

ft/sec 

34 

V(2,2,2)*12. 

A22 

in/ sec 

15 

PR 

p 

rad/sec 

35 

V(2.3,2)*12. 

A22 

in/ sec2 

16 

OR 

q 

rad/sec 

36 

V(2.1,3)*12. 

A23 

inches 

17 

RR 

R 

rad/sec 

37 

V(2,2,3)*12. 

A23 

in/ sec 

18 

PHI 

♦ 

rad 

38 

V(2,3,3)*12. 

A23 

in/ sec2 

19 

THETA 

0 

rad 

39 

FF(5)-FF(l)»ZCELL 

mh 

ft-lbs 

20 

AUXE(3) 

*1B 

non-dim 

40 

FF(4)+FF(2)*ZCELL 

lh 

ft-lbs 

21 


Figure  3.4  Comparison  of  P  with  dP/dt 


Comparison  ot  Q  with  dQ/dt 


Figure  3.6  Comparison  of  K  with  dR/dt 


Figure  3.7  Compa 


Rate  with  Urigina 


3.4.2  Data  Filtering 


The  REXOR  time  history  data  contain  high-frequency 
oscillatory  components,  with  a  frequency  of  four  per  revolution, 
exicted  by  the  rotation  of  the  four  blades.  This  frequency, 

4  x«R)  is  equal  to  26.94  Hz.  This  is  a  much  higher  frequency 
than  any  to  which  the  pilot  could  respond.  It  could  therefore  be 
classified  as  a  vibration  rather  than  a  response.  For  this 
reason,  it  was  decided  to  filter  out  these  higher  frequencies. 

To  achieve  this  filtering,  a  digital  implementation  of  a 
fourth-order  Butterworth  filter,  with  a  cutoff  frequency  of 
13.47  Hz  (2  x  was  used.  By  passing  the  data  through  the 

filter  twice,  once  forward  and  once  backward,  it  was  possible  to 
achieve  virtually  zero  phase  shift  throughout  the  pass  band. 

This  also  effectively  doubled  the  order  of  the  filter,  thus 
producing  a  steep  cutoff  above  the  filter  frequency.  Samples  of 
the  filtered  REXOR  data  for  the  hover  flight  condition,  main 
rotor  collective  blowing  input,  are  shown  in  Figures  3.9  through 
3.13. 
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Figure  3.1 
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0  Effect  of  Fourth-Order  Butterworth  Filter  ( “g* 
on  Vertical  Velocity  (Hover,  Collective  Input) 


IV.  IDENTIFICATION  OF  REDUCED-ORDER  MODEL 


4.1  OVERVIEW 

This  section  describes  a  methodology  for  extracting 
reduced-order  dynamic  models  of  the  X-wing  configuration  from 
REXOR  simulation  data.  The  data  processing  methodology  is  based 
on  an  adaptatin  of  system  identification  techniques  which  have 
been  used  to  process  flight  test  data  L  4 ] .  Even  though  the 
described  model  extraction  technique  is  demonstrated  with  the 
Lockheed  X-wing  configuration,  it  is  intended  to  be  applicable  to 
any  modeling  task  that  requires  the  extraction  of  reduced-order 
models  from  complex  simulations. 

The  scope  of  this  chapter  includes  a  discussion  of  system 
identification  as  it  applies  to  the  extraction  of  reduced-order 
models,  a  discussion  of  the  X-wing  modeling  approach,  and  a 
presentation  of  X-wing  reduced-order  modeling  results. 

4.2  SYSTEM  IDENTIFICATION  OVERVIEW 

System  identification  is  a  data  processing  technique  which 
has  been  used  to  estimate  the  parameters  of  a  math  model  from 
time  histories  of  the  dynamic  response  of  that  system.  System 
identification,  also  called  "parameter  identification,"  has  been 
used  for  a  variety  of  different  applications  (see  Table  4.1). 
Because  of  this  diverse  base  of  applications,  several  methods  of 
approach  for  identifying  moel  structures  and  parameter  values 
have  been  developed.  The  "best"  identification  algorithm/method 
is  the  one  which  is  best  suited  to  the  particular  requirements  of 
each  application. 

The  optimal  subset  regression  (OSR)  technique  is  well  suited 
for  the  task  of  extracting  reduced-order  models  from  the  REXOR 
simulation  generated  time  histories.  This  choice  is  based  on 


Role  of  Parameter  Identification  for  Aircraft 
Design,  Test,  and  Evaluation 

REDUCED-ORDER  MODEL  EXTRACTION  FROM  COMPLEX  SYSTEM  MATH  MODEL 
(e.g. ,  ROTOR  SIMULATION  SUCH  AS  REXOR) 

-  ADEQUATE  FIDELITY  MODEL  FOR  REAL-TIME  SIMULATION 
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FLIGHT  TEST  PLANNING 
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-  DATA  PROCESSING  SOFTWARE  DEVELOPMENT  AND  VALIDATION 
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AERODYNAMIC  MODEL  DOCUMENTATION  FROM  FLIGHT  TEST  MEASUREMENTS 

-  DATA  FORMAT  COMPATIBLE  WITH  USER  REQUIREMENTS 

-  PREFLIGHT  PREDICTION/FLIGHT  OBSERVATION  CORRELATION 

-  SIMULATION  MATH  MODEL  VALIDATION 

STATE  MODEL  EXTRACTION  FROM  FLIGHT  TEST  MEASUREMENTS 

-  APPLICATIONS:  AIRFRAME,  CONTROL  SYSTEM,  PROPULSION  SYSTEM,  ROTOR 
SYSTEM,  STRUCTURAL  MODES 

-  CONTROL  SYSTEM  PERFORMANCE  ANALYSIS  AND  VALIDATION 

-  HANDLING  QUALITY  PARAMETERS  FROM  IDENTIFIED  STATE  MODELS 


three  considerations.  First,  OSR  formulates  the  identification 
of  parameters  as  an  equation  error  problem.  With  the  equation 
error  technique,  the  parameters  are  the  coefficients  of  the 
independent  variables  that,  when  combined  in  a  linear  equation, 
form  the  dependent  variable  (i.e.,  P  =  KQ  +  K^p  +  K2Q  + 

K3A1B  +  ...).  Thus,  OSR  can  be  used  to  define  the  most 
significant  model  parameters  (i.e.  the  model  structure)  for  each 
degree  of  freedom.  The  aspect  of  obtaining  the  statistically 
most  significant  model  structure  is  discussed  in  more  detail 
later.  The  second  reason  for  selecting  the  OSR  algorithm  is  that 
it  is  computationally  a  very  efficient  program  to  use.  The  third 
consideration  in  the  selection  of  OSR  pertains  to  the  source  of 
the  time  history  data  which  is  used  to  identify  the  reduced -order 
model.  For  this  study,  the  data  were  generated  from  the  REXOR 
simulation  of  the  X-wing  configuration.  Theoretically,  if  the 
simulation  is  set  up  properly,  the  resulting  measurement  time 
histories  should  be  free  of  instrumentation  errors  and  noise. 

The  OSR  algorithm  works  best  with  uncorrupted  measurements.  As 
an  aside,  when  dealing  with  flight  test  measurements,  which  are 
often  corrupted,  a  more  sophisticated  identification  algorithm, 
such  as  the  maximum  likelihood  approach,  must  be  used. 

The  computational  steps  which  art  followed  by  the  OSR 
algorithm  are  described  as  follows. 

OSR  is  an  algorithm  which  adds  and  deletes  variables  to  a 
particular  model  in  an  iterative  manner.  The  algorithm  uses 
statistical  hypotehsis  testing  techniques  based  on  the  Fisher  F 
ratio  (e.g.  F-tests).  Formally,  this  ratio  measures  the 
difference  in  fit  error  with  the  current  model  relative  to  the 
error  due  to  noise  and  model  uncertainties.  An  "equation" 
F-ratio  measures  the  entire  model  fit  relative  to  error  and  a 
"parameter"  F-ratio  measures  the  incremental  improvement  in  fit 
due  to  addition  or  deletion  of  a  parameter  in  the  model. 

Starting  with  a  list  of  possible  variables,  the  algorithm  enters 
the  first  variable  with  the  highest  partial  correlation  to  the 


observations  (i.e.  the  dependent  variable).  The  contribution  of 
this  variable  to  reducing  the  fit  error  is  made,  and  a  new 
variable  is  entered.  Subsequent  tests  add  and  delete  variables 
to  improve  the  "fit."  This  procedure  is  repeated  until  no 
variables  are  entered  (or  removed)  from  the  model. 

The  model  structure  selected  from  the  OSR  results  is  based 
on  a  consideration  of  two  statistical  figures  of  merit  and  a 
comparison  of  candidate  model  structures  with  preflight 
predictions.  This  comparison  also  provides  a  means  to  validate 
the  selected  model  structure.  The  two  statistical  figures  of 
merit  are  the  multitude  correlation  coefficient  (R  )  and  the 
equation  F  ratio.  For  a  perfect  fit,  R  =  1.  Generally, 

R2  -+•  1  as  additional  terms  are  added  to  the  model.  The 
equation  F  ratio  relates  fit  goodness  to  fit  error  weighted  by 
the  degrees  of  freedom  of  the  model.  The  desired  model  is  one 
that  has  a  good  fit  (R2  1J  with  the  equation  F  ratio 
maximized. 

One  of  the  issues  addressed  at  the  outset  of  this  study  was 
the  feasibility  of  using  system  identification  to  define  the 
unaugmented  math  model  parameters  of  an  aircraft  employing  a 
high-gain  feedback  control  system.  The  question  arises  from  the 
fact  that  high-gain  feedback  loops  prevent  excitation  of  the 
unaugmented  aircraft  dynamics.  Previous  studies  have  shown  that 
the  unaugmented  dynamics  can  be  identified  from  the  augmented 
dynamics,  because  the  modes  suppressed  by  the  control  system  are 
reflected  in  the  control  system  output  (i.e.  the  control 
commands)  time  histories. 

Thus,  the  math  model  for  the  unaugmented  vehicls  is  obtained 
from  measurements  of  the  airplane's  response  and  the  control 
inputs.  For  the  X-wing  configuration,  the  control  inputs  are  the 
hub  moments,  collective  and  cyclic  blowing,  and  tail  rotor 


4.3  X-WING  MODELING  APPROACH 


Potential  applications  for  a  reduced-order  model  include  an 
assessment  of  handling  qualities,  synthesis  of  flight  control 
system  control  laws,  and  analysis  of  dynamic  coupling  between  the 
rotor/flight  control/propulsion  systems.  Linear  models  with  good 
physical  insight  are  desireable. 

The  desired  modeling  approach  is  based  both  on  the  intended 
useage  of  the  reduced-order  model  and  the  X-wing  configurational 
details . 

A  key  feature  of  the  X-wing  configuration  is  its  very  rigid 
rotor  system.  Because  of  the  rotor  blade's  structural  stiffness, 
most  maneuvers  produce  only  small  amounts  of  blade  deflection. 
Thus,  the  X-wing's  rotor  system  is  not  used  for  thrust  vectoring 
through  blade  flapping.  In  addition,  the  first  flapping  degree 
of  freedom  occurs  at  a  frequency  greater  than  the  principal  rotor 
frequency. 

Because  of  its  extreme  rigidity,  the  rotor  exerts  strong 
precessional  moments  on  the  aircraft  as  the  result  of  commanded 
angular  rates.  Coupling  moments  are  generated  by  aircraft 
rotations  orthogonal  to  the  rotor  rotation.  An  established  roll 
rate  produces  a  pitch  acceleration,  and  an  established  pitch  rate 
generates  a  roll  acceleration.  These  coupling  moments  are  thus 
primarily  dependent  on  rotor  moment  of  inertia  and  RPM. 

Additional  coupling  effects  due  to  rotor  blade  aerodynamics  arise 
from  nonuniformities  in  inflow,  downwash,  and  blade  loading. 

These  variables  are  a  function  of  rotor  RPM,  collective  and 
cyclic  blowing,  and  aircraft  airspeed  and  angle  of  attack. 
Fuselage  aerodynamic  characteristics  are  dependent  on  rotor  RPM, 
collective  and  cyclic  blowing  levels,  aircraft  angle  of  attack, 
sideslip  and  airspeed,  aircraft  rotational  rates  and  tail  rotor 
collective . 

A  linear  state  model  representation  is  the  most  appropriate 
form  for  the  reduced-order  model.  The  state  model  formulation  is 


ideally  suited  for  control  system  synthesis,  coupled  system 
analysis,  and  vehicle  handling  quality  analysis  problems.  The 
state  model  formulation  is  based  on  constant  coefficient 
equations  that  define  the  time  rate  of  change  of  the  state 
variables  as  a  function  of  current  values  of  the  state  and 
control  variables.  In  addition,  the  state  and  control  variables 
can  be  combined  algebraically  to  form  output  variables.  The 
equations  below  illustrate  the  state  model  formulation: 

x  =  Fx  +  Gu 

y  =  Hx  +  Du 

where  F,  G,  H,  D  are  constant  coefficient  system  matrices,  x 
is  the  state  vector  (e.g.,  axial  velocity,  body  pitch  rate,  rotor 
collective  flapping  position),  u-  is  the  control  vector  (e.g. 
collective  blowing,  cyclic  blowing,  etc.),  and  y  is  the  output 
vector  (e.g.  rotor  hub  pitch  torque).  As  a  side  comment  at  this 
point  in  the  discussion,  the  identification  of  the  elements  that 
comprise  the  system  matrices  (F,  G,  H,  D)  from  the  REXOR 
simulation  data  is  the  technical  objective  of  the  study 
documented  in  this  report.  Continuting  with  the  state  model 
formulation  discussion,  the  model  form  shown  above  can  be 
constructed  from  a  number  of  subsystem  state  models.  This 
concept  is  illustrated  in  Figure  4.1,  which  shows  the  make-up  of 
a  coupled  system  state  model  that  includes  the  body,  stability 
augmentation  system,  rotor,  drive  train,  engine  and  fuel 
controller  degrees  of  freedom.  The  scope  of  this  study  pertains 
to  only  the  airframe  and  rotor  degrees  of  freedom.  The  REXOR 
simulation  was  run  at  constant  rotor  RPM;  thus,  decoupling  the 
rotor  from  the  engine  and  transmission. 

The  X-wing  fuselage  plus  rotor  state  model  can  be  define  in 
one  of  two  possible  forms:  viz,  full-order  state  model  and  a 
quasi-static  representation.  For  the  full-order  formulation,  the 
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rotor  degrees  of  freedom  are  modeled  explicitly  (see  Figure 
4.2).  For  the  illustrated  example,  the  rotor  degrees  of  freedom 
are  described  by  the  tip  path  plane  displacements  and 
velocities.  For  the  quasi -steady  formulation,  the  rotor  states 
are  implicitly  modeled  by  fuselage  states.  The  procedure 
followed  in  reducing  a  full-order  state  model  to  a  quasi-static 
approximation  is  shown  in  Figure  4.3. 

The  quasi-static  formulation  was  selected  for  the  X-wing 
study  due  to  the  high  stiffness  of  the  rotor  blades.  Because  of 
this  stiffness,  the  flapping  modes  are  well  separated  from  the 
fuselage  degrees  of  freedom.  From  a  flight  control  engineer's 
point  of  view,  the  X-wing  rotor  behaves  like  a  high  bandwidth 
control  actuator. 

Specific  modeling  features  for  the  quasi-steady  X-wing  model 
are  described  as  follows: 

(1)  The  state  and  control  variables  are  referenced  to 
their  initial  values  (i.e.  X  =  %iQ  +  Fax  +  Gau). 

This  was  done  since  the  REXOR  simulation  was  not 
always  initialized  to  zero. 

(2)  Unsteady  aerodynamic  contributions  for  the  moment 
equations  were  modeled  in  terms  of  A  u,  av,  aw 
derivatives . 

(3)  The  angular  acceleration  equations  model  the  rotor 

by  a  hub  moment  plus  a  moment  resulting  from  the 
transfer  of  the  hub  force  to  the  center  of 
gravity.  The  hub  force  terms  are  proportional  to 
Ay  for  the  P  equation,  Ax  and  Az  for 

the  Q  equation  and  Ay  for  the  R  equation. 

(4)  The  inertial  coupling  between  P  and  R  due  to 
IXz  f  0  is  not  modeled  explicitly.  The  P 

•  • 

equation  does  not  include  R  terms,  and  the  R 

equation  does  not  include  P  terms. 

The  equations  of  motion  used  for  the  X-wing  analysis  are 
presented  in  Table  4.2.  The  equations  are  developed  for  a  body 
axis  system  representation,  modeled  the  kinematic  relationship 
with  perturbation  motion  and  assume  untrimmed,  wings-level 
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Figure  4.2  General  Rotorcraft  State  Model 


Table  4.2 

Equations  of  Motion  Used  for  X-Aing  Analys 


u  =  Ujc  +  R0Av  -  QqAw  -  wQ  aQ  +  v0AR 

-  g  cos  0Q  A0  +  Ax 

v  =  Vjq  -  Rqau  +  Pqaw  +  wqaP  -  uqaR 
+  g  cose  A8  +  a 

O  y 

W  =  +  QqAU  -  P0Av  -  V0AP  +  uqAQ 

-  g  sin  e  Ae  +  A 

o  z 


Q  =  QlC  + 


R 


<P  -  P  +  (Q  s in  ^>  +  R  cos  )  tan  e 
0  =  Q  cos  <P  -  R  sin  <t> 


initial  conditions.  The  applied  acceleratons  (i.e.  Ax  ,  Ay , 

K,  aMx/Ix,  AMy/Iy,  AMZ/IZ)  are  described  by  the 
identified  X-wing  math  model.  The  a  priori  model  structure, 
which  includes  all  possible  terms  for  the  model  structure 
determination  and  parameter  identification  analysis,  is  presented 
in  Table  4.3. 


4.  4  IDENTIFIED  MODEL  RESULTS 


Table  4.4  presents  a  summary  of  the  equation  fit  statistics 

(R2  and  Fgq)  for  the  identified  models  for  each  flight 

condition.  These  statistics  are  based  on  the  model  fit  errors 

(i.e.  the  difference  between  the  identified  model  and  the 

equation  value  calculated  from  REXOR  data)  for  all  four  test 

inputs  for  each  flight  test  condition.  R  ,  the  square  of  the 

equation  multiple  correlation  coefficient,  approaches  one  for 

reducing  fit  error.  F£q,  the  equation  F  ratio,  is  related  to 

R2,  and  it  places  a  penalty  on  the  number  of  terms  included  in 

2 

the  model.  The  relationship  between  R  and  Fgq  is  given 
by 


R2  N 
EQ  Y^R1  U 


1) 


where  N  =  number  of  data  points  modeled,  and 
M  =  number  of  parameters  in  the  model. 

The  identified  model  results  are  presented  in  Tables  4. 5 
through  4.12.  Each  of  the  tables  shows  the  regression  results 
for  one  dependent  variable  for  each  of  the  two  flight 
conditions.  Graphs  showing  a  comparison  of  the  identified  model 
and  the  original  time  history  are  shown  in  Figures  4.4  through 
4.35.  The  figures  show:  the  four  maneuvers  for  the  hover  flight 
condition,  an  enlarged  view  of  one  maneuver  for  the  hover  flight 
condition,  the  four  maneuvers  for  the  100-fps  flight  condition, 
and  an  enlarged  view  of  one  maneuver  for  the  100  fps  flight 
condition.  The  model  is  shown  as  a  solid  line,  and  the  original 
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Table  4.3 

X-Wing  A  Priori  Applied  Acceleration  Models 


Axial  Acceleration  (Ax) 

Ax  =  AXIC  +  Xu  Au  +  X„  Av  +  Xw  Aw  +  XD  AP  +  Xn  AQ  +  XD  aR 


'u  UV  AW  AP  ur  "  AQ  ^  ^  ar 


+  XA1B  M1B  *  XB1B  4B1B  *  XRp0  iRP0  *  Xe0TR  a80TR 


Lateral  Acceleration  (Ay) 

Ay  =  Ay jq  +  au  +  Yu  av  +  Yw  aw  +  Yp  aP  +  Yq  aQ  +  YR  aR 

+  YA1B  aA1B  +  YB1B  aB1B  +  YRp0  aRP0  +  Y0otr  a90TR 

Vertical  Acceleration  (A^) 

AZ  =  AZIC  +  Au  +  Zu  Av  +  Aw  +  ZP  aR  +  Zq  aQ  +  ZR  aR 

*  ZA1B  **1B  *  ZB1B  ‘SIB  +  ZRpo  8RP0  +  Z,0TR  A0OTR 

Roll  Acceleration  (Pi 


—j—  *  Lu  Au  +  Ly  AV  +  Lw  AW  +  Lp  AP  +  Lp  aQ  +  LR  AR 

X 

+  LA1B  AA1B  +  LB1B  AB1B  +  LRpo  ARP0  +  L0otr  a0otr 

+  LLH  a  LH  +  LMH  a  MH  +  Lu  A “  +  Lv  A  v  +  Lw  A  w  +  LAY  AY 
Pitch  Acceleration  (Q) 


—  =  Mu  AU  +  My  AV  +  Mw  AW  +  Mp  aP  +  Mp  aQ  +  MR  aR 


+  %  aA1B  +  M31B  aB1B  +  MRpo  aRP0  +  M90TR  a60TR 
+  MLH  aLH  +  MMH  aMh  +  Mj  AU  +  VI.  AV  +  M.  AW  +  MAX  Ax 


+  MAZ  AZ 


4  5 


Table  4.3  (Continued) 


Yaw  Acceleration  (R) 


~J~  =  Nu  Au  +  Nv  AV  +  Nw  AW  +  Np  aP  +  N0  aQ  +  NR  aR 


+  Na  AA,r  +  Nn  AB,r  +  Np  ARon  +  N„  A6nTD 
A1B  18  B1B  18  RP0  P0  90TR  0TR 


+  Nlh  ALH  +  Nmh  AMh  +  N*  Au  +  Av  +  N;  Aw  +  NAy  Ay 
Rotor  Hub  Rolling  Moment  (L^) 

lh  *  lhic  +  lhu  Au  +  lhv  Av  +  lhw  Aw  +  LHp  AP  +  LHq  aQ 

+LHrAR+LHa  aaib  +  lhr  abib+lhr  arpo 

aib  bib  rpo 

A9DTR  +  LH.  Au  +  Lu  Av  +  L„  Aw 
%TR  0TR  Hu  Hv  Hw 

MH  *  MHIC  +  MH  Au  +  MH  Av  +  MH  Aw  +  MH  aP  +  MH  aQ 

+  m  4r  +  n  aa18*m  4b1b+«  4rp0 


MH9qtr  490TR  +  \  A“  +  MH‘  4V  *  HH-  A" 
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Table  4.4 

Summary  of  Identified  Model  Equation  Error  Fit  Statistics 


EQUATION 


HOVER 

100 

fps 

R2 

feq 

R2 

feq 

.998 

96,555 

.905 

1,840 

.999 

181,550 

.973 

6,326 

,999 

506,990 

.999 

541,410 

.904 

1,207 

.951 

2,686 

.919 

1,451 

.977 

4,459 

.996 

31,106 

.997 

43,888 

.994 

22,666 

.989 

13,394 

.997 

50,886 

.997 

67,052 

f 


time  history  is  shown  as  a  series  of  circles.  Not  all  points 
used  were  plotted,  so  some  of  the  apparent  mismatch  during 
high-frequency  transients  is  due  to  the  plotting. 

Tne  data  presented  in  Table  4.4  indicate  that  the  identified 
reduced -o rde r  models  reproduce  the  applied  accelerations 
generated  by  REXOR  for  each  of  the  four  control  inputs.  The  fit 
errors  are  very  good  (i.e.  >  .99)  with  the  exception  of  P 

and  Q  at  both  flight  conditions,  and  Ax  and  Ay  at  the 
100-fps  flight  condition. 

The  time  histories  for  Ax  and  Ay  (see  Figures  4.6 
and  4.10)  indicate  that  the  modeling  problem  with  these 
accelerations  is  primarily  with  the  high-frequency  portion  of 
collective  blowing  input.  The  models  did  not  account  for  load 
sharing  between  the  main  and  tail  systems,  and  it  also  did  not 
include  direct  thrust  terms. 

Examination  of  the  regression  statistics  and  graphical 
comparisons  showed  that  a  modeling  discrepancy  existed  in  the 
pitch  and  roll  accelertions  and,  to  a  lesser  extent,  the  pitch 
hub  moment.  The  angular  acceleration  errors  may  be  related  to 
the  discrepancy,  previously  noted,  between  the  pitch  and  roll 
rates  and  accelerations  from  REXOR.  Comparison  of  the  model 
error  for  Q  at  hover  (Figure  4.21)  with  the  original  data  for 
at  hover  (Figure  4.33),  both  with  the  same  input 
(longitudinal  cyclic),  shows  a  strong  similarity.  These 
comparisons  indicate  that  the  pitch  acceleration  (Q)  and  the  hub 
pitching  moment  from  REXOR  are  not  completely  consistent  with 
each  other. 
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Table  4 .  5 
Regression  on 


_  IDENTIFIED  MODEL 

oooooooooo  REXOR  SIMULATION  DATA 


IDENTIFIED  MODEL 

oooooooooo  REXOR  SIMULATION  DATA 


SECONDS  I 


_  IDENTIFIED  MODEL 

oooooooo oo  REXOR  SIMULATION  DATA 


SECONDS)  |  TIME  (SECONDS) 


IDENTIFIED  MODEL 

oooooooooo  REXOR  SIMULATION  DATA 


Forward  Acceleration  (100  fps.  Collective  Input) 


Tab  Le  4.6 
Regression  on 


IDENTIFIED  HOOEL 


TIME  (SECONDS)  |  TIME  (SECONDS) 
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oooooooooo  REXOR  SIMULATION  OATA 


Figure  4.9  Lateral  Acceleration  (Hover,  Collective  Input) 
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Table  4.7 
Regression  on 


_  IDENTIFIED  MODEL 
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Figure  4.12  Normal  Acceleration  (Hover) 
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Figure  4.13  Normal  Acceleration  (Hover,  Collective  Input) 
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Table  4.8 
Regression  on  P 


HOVER 

100  FPS 

VARIATION  EXPLAINED  8Y  R2: 

0.90392160 

VARIATION  EXPLAINED  BY  R2: 

0.95128113 

EQUATION  F-RATIO:  1207.4 

EQUATION  F-RATIO:  2686.2 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

VARIABLE 

COEFFICIENT 

PARAH.  F-RATIO 

Au 

0.0207611 

64.75 

-0.0611522 

225.7 

Av 

0.0102289 

103.3 

0.0301222 

331.9 

AW 

0.0206216 

7.077 

-0.0948021 

60.42 

AP 

8.97866 

660.2 

aP 

-11.0907 

1310. 

AQ 

47.0976 

10290. 

aQ 

70.1132 

5325. 

AR 

-1.88311 

4523. 

iAlB 

4.21199 

334.0 

CO 

H 

CO 

<1 

1.41577 

115.0 

iBlB 

1.66284 

314.2 

iRP0 

0.381965 

15.59 

ARpO 

-5.75918 

24.07 

a0OTR 

7.35574 

367.3 

ie0TR 

-1.08932 

9.330 

alh 

0.00424158 

14560. 

alh 

0.00454312 

14220. 

-0.000251747 

107.7 

AMh 

-0.000182272 

90. 

a'u 

0.0115351 

29.94 

aw 

-0.068116 

121.2 

a'v 

-0.00755679 

73.71 

ay 

1.90513 

5794. 

ay 

1.10465 

2112. 

CONST. 

-1.30088 

1914. 

CONST. 

-2.19794 

1096. 

1 

I 

i 
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Figure  4.16  Roll  Acceleration  (Hover) 
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Table  4.9 
Regression  on  Q 


HOVER 

100  FPS 

VARIATION  EXPLAINED  BY  R2: 

0.91874158 

VARIATION  EXPLAINED  BY  R2- 

0.97661561 

EQUATION  F-RATIO:  1451.0 

EQUATION  F-RATIO:  4459.4 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

1* 

0.0375316 

1093. 

0.0375001 

178.6 

AP 

-12.6538 

5434. 

0.0146625 

159.3 

AQ 

5.26030 

716.1 

Aw 

0.0543127 

31.78 

AR 

-0.287913 

388.1 

AP 

-14.1868 

3622. 

AA1B 

-2.42522 

1971. 

AQ 

-14.2434 

26.49 

A81B 

-4.66679 

3069. 

AR 

14.3501 

244,5 

iRP0 

-1.14406 

124.0 

AA1B 

-0.560806 

11.04 

A9  OTR 

1.53243 

127.1 

Abib 

-0.933818 

377.3 

A4 

-0.000144149 

65.89 

ARro 

-5.26926 

699.6 

amh 

0.00103564 

5466. 

i60TR 

2.00070 

96.73 

4  u 

0.0294313 

527.1 

ALH 

0.000372978 

507.4 

A  v 

-0.0151663 

868.3 

amh 

0.000873825 

4357. 

Aw 

0.0245589 

267.6 

A  li 

0.0224886 

22.02 

AX 

1.29859 

3415. 

A  v 

0.147567 

237.7 

CONST. 

-0.0302397 

25.81 

A  w 

0.179023 

41.30 

Ax 

1.41308 

980.5 

Az 

-0.124869 

19.94 

CONST. 

-2.89699 

_ 

9.815 
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Figure  4.21  Pitch  Acceleration  (Hover 
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Table  4.10 
Regression  on  R 


HOVER 

VARIATION  EXPLAINED  BY  S2: 

EQUATION  F-RATIO:  31106. 

0.99614894 

VARIATION 

EQUATION 

100  FPS 

EXPLAINED  8Y  R2: 

-RATIO:  43888. 

0.99708450 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

Au 

-0.00344913 

10.1 

Au 

0.00858504 

97.34 

Aw 

0.0220235 

54.33 

Av 

-0.0132009 

1890. 

AP 

-1.59912 

159.9 

Aw 

-0.0530431 

476.3 

AQ 

-9.30556 

3232. 

AP 

2.05455 

1273. 

AR 

0.544170 

2249. 

AQ 

-12.0759 

4156. 

“A1B 

3.74247 

5521. 

ar 

1.51370 

15.37 

i8lB 

-1.37695 

1412. 

*aib 

-0.201839 

11.46 

ARpO 

3.97052 

4284. 

iSlB 

-0.166797 

30.39 

i90TR 

-2.49114 

274.0 

ARpo 

-0.673275 

90.53 

iLH 

-0.000863718 

4588. 

-0.000792314 

11570. 

ahh 

0.000168149 

373.6 

amh 

0.0000232729 

39.57 

A  U 

0.00697193 

68.90 

■ 

0.0166052 

17.50 

A  v 

0.00975561 

781.5 

gpiP 

0.00948151 

63.59 

A  w 

-0.0156250 

294.5 

-1.52544 

83330. 

Ay 

CONST. 

-1.35129 

3.77182 

20380. 

20700. 

1.03957 

25040. 

*  1 


CONST 


3.77182 


20700 
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Table  4.11 
Regression  on 


HOVER 

VARIATION  EXPLAINED  3Y  R2: 

EQUATION  F-RATIO:  22666. 

0.99350286 

100  FPS 

VARIATION  EXPLAINED  BY  R2 

EQUATION  F-RATIO:  13394. 

:  0.98905438 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

AU 

6.26771 

15.20 

AU 

8.89394 

11.87 

AV 

-24.9992 

800.7 

AV 

-20.8489 

493.3 

AW 

-11.6902 

.  7.785 

AW 

40.5570 

60.31 

AP 

-1262.32 

41.95 

AP 

-350.810 

3.248 

AQ 

-5160.66 

448.2 

AQ 

-16525.1 

1774. 

AR 

563.580 

2777. 

aa1B 

226.245 

26.56 

iAlB 

3112.06 

699.7 

ARpO 

4675.30 

563.8 

A81b 

997.119 

229.7 

a90TR 

5229.36 

254.8 

ARpo 

586.745 

126.3 

A  u 

45.5093 

61.29 

49otr 

2200.32 

87.50 

a; 

-5.48350 

145.4 

Au 

-14.3104 

125.3 

Aw 

39.3007 

885.0 

Ay 

-417.352 

1203. 

ay 

-599.450 

1397. 

CONST. 

260.987 

55.03 

CONST. 

-722.783 

2404. 
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Figure  4.28  Roll  Hub  Moment  (Hover) 
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Figure  4.29  Roll  Hub  Moment  (Hover,  Collective  In 
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Figure  4.31  Roll  Hub  Moment  (100  ips,  Collective 


Table  4.12 
Regress  ion  on  MH 


HOVER 

VARIATION  EXPLAINED  BY  R2: 

EQUATION  F-RATIO:  50886 

0.99748439 

100  FPS 

VARIATION  EXPLAINED  BY  R2: 

EQUATION  F-RATIO:  67052 

0.39712998 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

VARIABLE 

COEFFICIENT 

PARAM.  F-RATIO 

Au 

24.9415 

40.95 

Au 

37.4316 

99.46 

AV 

-36.5539 

553.6 

A  v 

-14.4361 

147.6 

Aw 

-23.6087 

9.556 

AP 

-7307.50 

470.6 

AP 

-3263.56 

186.3 

AQ 

5160.67 

82.55 

AQ 

4914.21 

148.1 

"S1B 

-52.6170 

13.31 

AR 

-78.2995 

22.30 

ARP0 

-3103.93 

176.2 

iAlB 

4368.84 

251.8 

iSQTR 

2592.14 

245.8 

a8ib 

1425.29 

108.3 

Aw 

-29.6058 

52.11 

ARpo 

1924.18 

157.4 

AX 

245.398 

46.03 

A9OTR 

1255.79 

82.09 

CONST. 

2088.23 

9162. 

AU 

9.51980 

21.90 

Aw 

108.436 

84.95 

AX 

2594.15 

2426. 

Az 

-112.266 

46.03 

CONST. 

-2726.96 

49.35 
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igure  4.33  Pitch  Hub  Moment  (Hover,  Longitudinal  Cyclic  Input) 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

In  general,  the  identified  models  of  the  X-wing 
configuration  show  good  correlation  with  the  original  data. 
Further,  the  identified  models  are  all  simple  and  of  much  lower 
order  than  those  used  in  the  REXOR  simulation.  This  indicates 
that  it  is  indeed  feasible  to  use  system  identification 
techniques  as  a  method  of  developing  a  low-order  math  model  for  a 
rotorcraft  from  a  simulation  math  model  that  represents  the  rotor 
by  its  individual  blades  in  a  rotating  frame  of  reference. 

The  discrepancies  seen  in  the  identified  results  seem  to  be 
related  to  the  discrepancies  noted  in  the  data  produced  by  the 
REXOR  program  (see  Sections  3.3  and  3.4  of  this  report),  although 
time  did  not  permit  a  thorough  analysis  of  these  d i scrpanc i es . 
This  should  serve  as  a  caveat  in  the  use  of  the  identified  model, 
as  well  as  the  REXOR  X-wing  simulation  as  implemented  at 
NASA -Ames . 

5.2  RECOMMENDATIONS 

In  order  to  fully  treat  the  validity  of  the  identified 
model,  a  stand-alone  simulation  using  this  model  should  be 
compared  to  the  filtered  REXOR  time  histories.  In  addition,  the 
REXOR  model  itself  should  be  verified,  and  the  causative  factors 
in  the  noted  data  inconsistencies  should  be  determined  and 
corrected.  This  was  not  done  due  to  time  constraints  and  the 
fact  that  no  documentation  exists  for  the  X-wing  version  of  the 
REXOR. 

In  general,  it  can  be  recommended  that  system  identification 
be  used  as  a  tool  for  developing  low-order  models  of  rotorcraft 
in  general.  This  may  be  done  from  simulation  data  or  from  flight 
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